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Graphene, being one-atom thick, is extremely sensitive to the presence of adsorbed atoms and molecules and, 
more generally, to defects such as vacancies, holes and/or substitutional dopants. This property, apart from 
being directly usable in molecular sensor devices, can also be employed to tune graphene electronic properties. 
Here we briefly review the basic features of atomic-scale defects that can be useful for material design. 
After a brief introduction on isolated Pz defects, we analyse the electronic structure of multiple defective 
graphene substrates, and show how to predict the presence of microscopically ordered magnetic structures. 
Subsequently, we analyse the more complicated situation where the electronic structure, as modified by the 
presence of some defects, affects chemical reactivity of the substrate towards adsorption (chemisorption) of 
atomic/molecular species, leading to preferential sticking on specific lattice positions. Then, we consider the 
reverse problem, that is how to use defects to engineer graphene electronic properties. In this context, we 
show that arranging defects to form honeycomb-shaped superlattices (what we may call "supergraphenes") 
a sizeable gap opens in the band structure and new Dirac cones are created right close to the gapped region. 
Similarly, we show that substitutional dopants such as group IIIA/VA elements may have gapped quasi- 
conical structures corresponding to massive Dirac carriers. All these possible structures might find important 
technological applications in the development of graphene-based logic transistors. 
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Graphene, thanks to its extraordinary electronic and 
mechanical properties, is a potential candidate for a num- 
ber of applications. Being one-atom thick, it is extremely 
sensitive to the presence of adsorbed atoms and molecules 
(either physisorbed or chemisorbed on the surface) and, 
more generally, to defects such as vacancies, holes and/or 
substitutional dopants. This property, apart from being 
directly usable in molecular sensor devices, can also be 
employed to tune graphene electronic properties. 
In this Chapter we review those basic features of atomic- 
scale defects that can be useful for material design. Af- 
ter a brief introduction (Section |ll| of the main prop- 
erties determining the peculiar electronic structure of 
graphene, and the experimental realisation of defective 
substrates (Section III), we focus in Section IV on iso- 
lated "p^ defects" such as atom vacancies or adsorbed 
species which covalently bind carbon atoms. In particu- 
lar, we discuss in detail the formation of so-called midgap 
states and the microscopically ordered magnetic struc- 
tures which give rise to. In Section [V| we analyse the elec- 
tronic structure of multiple defective graphene substrates 
and show, in particular, how it is possible to use simple 
rules to predict the presence of magnetic moments and 
midgap states by looking at the defect locations on the 
lattice. Subsequently, we analyse the more complicated 
situation where the electronic structure, as modified by 
the presence of some defects, affects chemical reactivity 
of the substrate towards adsorption (chemisorption) of 
atomic/molecular species, leading to a preferential stick- 
ing on specific lattice positions. In Section[VT|we consider 
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the reverse problem, that is how to use defects (vacan- 
cies, adsorbed species, substitutional dopants, etc..) to 
engineer graphene electronic properties. This is possible 
nowadays since recent advances in lithographic and self- 
assembling techniques allow one to produce well-ordered 
structures and thus 'tune' the electronic bands. In this 
context, we show for instance how it is possible to open a 
band-gap in graphene and preserve at the same time the 
pseudo-relativistic behaviour of its charge carriers. We 
further analyse the case of substitutional dopants (group 
IIIA/VA elements) which, if periodically arranged, may 
show a gapped quasi-conical structure corresponding to 
massive Dirac carriers. All these possible structures 
might find important technological applications in the 
development of novel graphene-based logic transistors. 



II. THE TT-ELECTRON GAS 

Carbon atoms in graphene are arranged to form a hon- 
eycomb lattice tightly held by strong a bonds between 
sp^ orbitals which form occupied a bands at energies well 
below the Fermi level. The remaining valence electrons 
(one for each carbon atom) populate a tt band which lo- 
calises above and below the lattice with a node on the 
surface plane. An 'antibonding' tt* band is empty when 
the system is at T = K and charge-neutral, but can 
easily be occupied, e.g. by applying a gate potential 
in a typical field-effect transistor (FET) configuration. 
Such tt/tt* band system governs the low-energy (say up 
to ^ 2 eV) behaviour of charge carriers in graphene and 
is responsible for most of the extraordinary properties of 
this material. This "tt cloud" is the focus of this section, 
where we introduce the main theoretical tools used in 
this Chapter. 
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where a] ^ a) creates an electron with spin a =t, i 
on the i-th lattice site of the A(B) sublattice, the first 
two sums run over nearest neighbouring sites (ti is the 
hopping energy) and the second ones over sites which 
are nearest neighbours in each sublattice (^2 is the cor- 
responding hopping) In absence of magnetic fields the 
hoppings can be chosen real, and the accepted value 
for ti is ^ 2.7 eV while 1^2 1 << ti depends on the 
parametrization used. Neglecting overlap between or- 
bitals on different C atoms, the usual anticommutation 
rules [cj ^, c'j,cr']+ = ^c,c'^ij^a,a' (c = &) hold; hence, 
introducing the Fourier transformed operators ak,cr ac- 
cording to 
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where the sum runs over k points in the first Brillouin 
zone (BZ) (analogously for 6k, cr) the above Hamiltonian 
can be rewritten as 



^^^ = -ti^/(k)4^^6k,a + /^.c. 
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A. Tight-binding Hamiltonian 



or, in matrix notation. 



t2g{k) hf{k) 
hf*(k) t2g{k) 
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In building up a simple, one-electron model for these tt 
electrons only one writes the one-electron wave function 
as a linear combination of two Wannier basis functions 
built with Pz orbitals, one for each sublattic^, 

V^k(r) = CAV^i^(r)+CBV^^(r) 



^k(r) = -^E^"^^'^^^(^-I^/) 

(and similarly for t/j^^r)) where the sum runs over lat- 
tice vectors Rj within a large supercell S including N 
graphene unit cells and = Hj -\- Sa is the position of 
A site in the j-th cell Equivalently, in second-quantized 



Here /(k) and ^(k) are 'structure factors' for the 
nearest- and next-nearest neighbours, 

/(k) = J2 e"'*^'- 

i=l,3 

5(k) = ^ e-^^< 
Diagonalization is trivial and gives the energy bands. 



e(k)± = -t2g{k)±h\f{k)\ = -t2g{k)±h^/3 + g{k) 

(2) 

where |/(k)p = 3 + ^(k) has been used and the minus 
(plus) sign solution correspond to the tt (tt*) band (see 
e.g. Ref.sHEEj). Close to the i^(i^O Point |/(K + q)p - 



FIG. 1. Left panel: graphene unit cell (ai,a2), along with the vectors joining nearest- and next-to nearest neighbours, Si and 
Si respectively. Also indicated the position vectors Sa and Sb for A and B sites. The Wigner-Seitz cell is one the hexagons. 
Right: first Brillouin zone with the highest symmetry point indicated. The arrows are the reciprocal lattice vectors (ai,a2). 



v^q'^ and the dispersion is conical, giving rise to the so- 
called Dirac cones. Here vp = ^tia = where d is 
the carbon-carbon distance, ~ 1.42 A, and a the lattice 
constant. Consequently, the density-of-states (DOS) is 
linearly vanishing at zero energy, p{e) ~ 4|e|/7rA/3^i (per 
cell, with spin and valley degeneracies included), one of 
the fingerprints of massless Dirac electrons. Its vanishing 
value challenges one's intuition since experiments find a 
finite, non-zero minimum conductivity at this energy^. 

Albeit simple, this tight-binding model is accurate 
enough to correctly represent graphene tt bands, at least 
close to the high symmetry points K and K' . The lat- 
ter control the low-energy physics of charge carriers, and 
are the source of the exceptional interest in graphene. 
If only nearest-neighbours interaction is allowed the two 
sublattices form two disjoint sets where A- type sites con- 
nect to B-type sites only and vice versa. The Hamiltoni- 
ans is said bipartitic and displays an interesting symme- 
try: for each non-zero energy level e and eigenfunction 
= lA) -\- cb \B) (where \A)/\B) is non-zero on 
A/B lattice sites only), there exists a 'conjugate' level 
with energy — e and wavefunction |?/;_) = ca\A) —cb \B) . 
This is called electron-hole (e — h) symmetry since at half- 
filling (as it is case of graphene with one electron per 
site), the Fermi level lies at zero energy, and the above 
symmetry relates electron and holes. For a proof, just 
apply a phase-change to one of the two sets of sublat- 
tices stated, e.g. bi^a —bi^a-, as this converts H into 
-H. 

Electron-hole symmetry, as we shall see in the following, 
plays an important role in graphene, even if it holds only 
approximately {i.e. with nearest-neighbour interaction 
only and neglecting orbital overlap). Here we just notice 
that, because of such symmetry, the band-structure is ex- 
pected to have a gap at the Fermi level unless there are 
specific reasons for having energy levels exactly at zero. 
As we now show, the specific reason is provided by the 
spatial symmetry of the substrate. 

Graphene lattice is highly symmetric. Its Wigner-Seitz 
unit cell has the same point symmetry of benzene, namely 
it belongs to the Dq^ point group, see figj2j which is 



the point group for symmetry operations in real-space. 
For Bloch electronic states with /c-vector k, symmetry 
is reduced to that subgroup of Dq^ which either leaves 
k invariant or transform it into one of its images, i.e. 
k ^ k + G with G a reciprocal lattice vectoil^. Such sub- 
group is known as /c-group at k, G(k), and determines the 
possible symmetry of the electronic states. For instance, 
at the K point the /c-group is D^h since only three-fold 
rotation axes and cf^ planes transform the K images into 
themselves^. A full analysis of the symmetry properties 
of Bloch electrons is given in figjsj left panel, where the k- 
groups are color-coded, grey for Cg, black for C2V1 green 
for D2hi blue for D^^h and red for DQh- The main point 
here is that graphene is sufficiently symmetric that allows 
/c-groups supporting ti(;o-dimensional irreducible repre- 
sentations {E irreps), namely Dq^ at F and D2,h at i^, K' . 
As spatial symmetry is (almost) compatible with e — h 
symmetry, a zero energy state results whenever the elec- 
tronic wavefunctions span a two-dimensional irreducible 
representation (odd in number, in general), i.e. they give 
rise to a doubly-degenerate level. This is exactly the case 
of the K {K') point, where Wannier functions built with 
Pz orbitals of the A and B sublattice span the E" irrep of 
the above D^^ A^-group . Notice also that this symmetry 
argument is enough to explain the conical dispersion of 
the energy at the K {K') point which makes graphene 
so attractive: without an inversion center, degeneracy is 
lifted already at first order in k • p perturbation theory 
when moving away from the BZ corners. 
While spatial symmetry is exact, e — h symmetry holds 
in the nearest-neighbour approximation only. Neverthe- 
less, since inclusion of higher order hopping terms does 
not modify the level ordering {i.e. the minimum of the 
TT* band lies always above the top of the tt band) the 
Fermi level at charge neutrality matches exactly the en- 
ergy where the E irrep is found (called Dirac point as it 
is the cone apex). 



4 




FIG. 2. Left panel: Point symmetry elements in graphene lattice. Reflection planes orthogonal to the page, cr^ and cr^, are 
replicated by the six-fold rotation axis Ce, along with the two-fold rotation axis on the page plane, C2. / is the inversion center, 
at the center of the Wigner-Seitz cell (solid line). The page plane is, of course, a reflection plane (cr^). Right panel: symmetry 
elements of the A:-group at the K point {Dsn). Black dots mark the K point and its images. 
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FIG. 3. Left panel: point-symmetry classification of Bloch states in graphene. /c-groups are color coded as indicated. Right 
panel: symmetry labels for graphene electronic states on the highest-symmetry lines of the BZ. 



B. Hubbard Hamiltonian 

The tight-binding (TB) Hamiltonian is a model in 
which each electron moves independently from the others. 
Despite it represents a good approximation for graphene 
energy spectrum, such a simple picture will necessarily 
fail in computing spin properties in all but the simplest 
situations. A simple way to include electron-electron in- 
teractions is given by the Hubbard model 

H = H^''^uY,<t<i (3) 

i 

where H^^ is the tight-binding Hamiltonian of eq{l| the 
sum runs over all carbon sites and fii^a- = a^i,o- are the 
corresponding number operators. This Hamiltonian com- 
bines the tendency of electrons to delocalize onto the lat- 
tice due to their kinetic (hopping) energy together with 
an "on-site" Coulomb repulsion that tends (for U > 0) to 
localize them to minimize double orbital occupation. 
The Hubbard model is a very useful tool for the study of 
magnetism in complex materials. It has long been used 
in the chemical communit}'^ -and proved to be rather 
accurate for such systems- to investigate excitation spec- 
tra in polycyclic aromatic hydrocarbons (today, graphene 
dots). Though simple, the model requires quite a large 



effort for its solution. Therefore, one often resorts to its 
mean-field approximation, 

H-^f ^JjTB^jjY^ ^.^^ ^^^^^^ ^jjY^ ^^.^^^ ^^^^ 

i i 

i 

where the average occupation number of one spin-species 
at a given site tunes an effective on-site energy for the 
other spin-species, e.g. e^^^ = +C/ {fii^i). This is essen- 
tially equivalent to an (unrestricted) Hartree-Fock ap- 
proach to the TT electrons and is useful, as compared 
with density-functional-theory (DFT) methods applied 
to the exact Hamiltonian, to study very large systems, of 
dimension comparable to those experimentally realized. 
Though we will not solve the Hubbard model in the fol- 
lowing, there are some exact, analytic results that can be 
obtained from it and that turn out to be important tools 
in discussing defects in graphene. 

C. Valence Bond picture 

An alternative, easy-to-use way of looking at graphene 
electronic structure is provided by the 'chemical picture'. 
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FIG. 4. The five possible perfect pairing Rumer diagrams for 
the benzene molecule with their correspondence with Kekule 
(1-2) and Dewar (3-5) resonance structures. 



With this we mean the traditional picture of chemical 
bonds as given by the Lewis structures and modified to 
account for the 'chemical resonance'. In this picture, elec- 
trons are mostly localized in atomic orbitals (usually hy- 
bridized) of the atoms forming the molecule, and couple 
in singlet pairs to form bonds and lone-pairs. For carbon 
atoms in graphene the three sp'^ orbitals (with one elec- 
tron each) are singlet-coupled with electrons in sp'^ or- 
bitals of neighbouring sites. The remaining electron (the 
one described by the TB Hamiltonian above) can couple 
with its counterpart of one of the three neighbours. The 
state of the system is a superposition of these different 
ways of binding, and the system gains energy from such 
a resonance phenomenon. 

This naive picture finds its root in the Valence Bond 
(VB) theory of chemical bond, which developed from 
the Heitler-London study of the H2 molecule, soon af- 
ter the foundation of quantum mechanics. The theory, 
as intensively pushed forward by Slater and Pauling, is 
a practical way of looking at the chemical bond and at 
the bond-breaking, bond-forming processes which are es- 
sential for chemical reactivity. It can also be turned 
into a variational method for the many-electron problem 
which uses a correlated wavefunction ansatz and captures 
the important part of the electron correlation^ In 
many respects, it has to be considered complementary to 
the Molecular Orbital (MO) approach, though the latter 
proved to be numerically more efficient. 
Valence Bond theory focuses on spin and builds the 
singlet wavefunction of an even-numbered ground-state 
molecule as a 'product' of singlet pairs, one for each bond 
(pairs of orbitals), thereby identifying a chemical for- 
mula. For less standard species such as graphene, differ- 
ent products are equally likely and the correct wavefunc- 
tion is the linear combinations of all the possible struc- 
tures. For instance, let us look at the benzene molecule 
as prototypical case of aromatic compounds. Consider- 
ing only the six tt electrons localized in their respective Pz 
orbitals, the possible linearly independent ("perfect pair- 
ing") functions can be schematically depicted as in figure 
|4j for six electron and an overall singlet state there are 
five couplings^^. In the graphical representation of Figjij 
known as Rumer diagrams, one represents each atomic 
center with a dot and uses a line for singlet coupling be- 



tween them. The chemical picture (bottom row of the 
same figure) uses only the 2 Kekule structures on the 
left, since the other 3 so-called Dewar structures gives a 
negligible contribution to the energy, as can be guessed 
from the bond pattern. 

Notice that VB theory is closely related to the Hubbard 
model discussed in the previous section: the atomic or- 
bitals housing the electrons need not be those of the free 
atoms. If they are 'polarized' by the environment {e.g. 
they are linear combinations of free-atomic orbitals) the 
VB ansatz accounts both for localization and band-like 
behaviour, as in the Hubbard model. It is not hard 
to show, indeed, that the Hubbard model for the H2 
molecule can be obtained from a simpl^^ VB ansatz to 
the two-electron wavefunction. In the following we will 
make a qualitative use of this chemical picture, as it pro- 
vides insights into defect formation and properties; see 
also Ref. [16| for its role in interpreting details of STM 
pictures. 

III. DEFECT FORMATION 

In this Section we describe defect formatiorJ^. As 
we shall see in the following, even though vacancies and 
adatoms turn out to have very similar effects on the tt 
electron system, we distinguish them here according to 
their preparation methods, i.e. high (e~, ions, etc.) vs. 
low (neutrals) energy beams. The reason is that only 
in the first case defects can be considered randomly ar- 
ranged. Adatoms at all but very low concentration tend 
to cluster on the surface, and understanding this phe- 
nomenon requires knowledge of how the electronic struc- 
ture of the substrate is modified upon formation of the 
very first defects. 

Notice, however, that defects are also naturally present 
in graphene as in any common material. Adatoms as hy- 
drogen or hydrocarbons, for instance, can be introduced 
by the preparation method, while other point defects af- 
fecting transport properties, e.g. charge-impurities, lo- 
cal potentials, etc., may result from imperfections on the 
substrates where graphene is accommodated. 

A. Electron and ion bombardment 

The irradiation by high energy particles is the main 
tool for creating defects in graphene and in other car- 
bon nanostructures. When the projectile particle im- 
pinges on the structure it transfers energy to the lattice. 
In bulk materials {e.g. in graphite) energy dissipation 
is rather effective, up to eventually stop the projectile, 
and it occurs through nuclear and electronic stopping 
mechanisms. Nuclear stopping is due to the collisions 
between the projectile and the carbon nuclei, an essen- 
tially classical process governed by momentum transfer 
and Coulomb interaction. On the other hand, (inelastic) 
electronic stopping occurs by the many possible electron 
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transitions in the material, hence promotion into conduc- 
tion band (hot electrons), ionizations, but also through 
plasmon excitations, photoemissions, etc. The relative 
importance of the two mechanisms depends on the beam 
energy, on the projectile mass and on the electronic struc- 
ture of the target material. Nevertheless a microscopic 
theory of energy dissipation in nanostructures is still un- 
der study since the models developed for bulk materials 
cannot be easily applied in a reduced dimensionality ma- 
terial such as graphene^^. 

The mechanism for the defect formation has been studied 
intensively in the last decade. In brief, when the energy 
transferred to an atom is larger than the so-called dis- 
placement threshold (~20 eV in case of graphite) this 
can leave its equilibrium position and move trough the 
bulk to form, for instance, a Frenkel pair or, for single 
layer graphene, a vacancy. Large ions can produce mul- 
tiple vacancies up to small holes in the lattice depending 
on their size. Electron beams produced in transmission 
electron microscopes (TEM) can instead be focused down 
to scales comparable to the carbon-carbon distance, giv- 
ing a precise control of the induced damage up to form 
single vacancies. Moreover TEMs allow a real-time imag- 
ing of the damage process and of t he ch emical reaction 
that follows the vacancies format ioiP^'^. 
The formation of a single vacancy in graphene leaves 
three a dangling bonds and it removes a tt electron. The 
first span a low-energy, one-dimensional irreducible rep- 
resentation [A irrep in the following) of the (local) D^^ 
point group and an E irrep. Therefore, the ground-state 
is degenerate and undergoes a Jahn- Teller distortion: the 
closure of two dangling bonds to form a pentagon, with 
an energy gain of about 0.2 eV. The strain induced by 
the other hexagons in the lattice prevents further distor- 
tions of the third unsaturated atom out of plane^^ and 
the final magnetic moment for such a structure has been 
reported to be between 1.0 and 1.5 /i^^^ localized on 
the unpaired site. When exposed to a hydrogen flux, 
the vacancy rapidly saturates its dangling bonds, with H 
atoms pointing slightly out of the graphene plane^. 
In the case of neutral-atom bombardment, the projectile 
can also react to form a covalent bond with a carbon 
atom. This is what happens by irradiating samples with 
low energies hydrogen atoms. It has been shown that 
at very low densities the chemisorbed H atom defect pro- 
duces STM images very similar to the single vacancy case. 
As already mentioned, at higher densities H atoms tend 
instead to cluster in dimers or larger structures due to 
electronic effects that will be discussed in the following. 
Nevertheless, when considering tt electrons only, vacancy 
and singly-bond chemisorbed species are equivalent, since 
a single electron is removed from the aromatic network 
of graphene. As an example, chemisorption of a single H 
atom is detailed in the following section. 



B. Sticking of atomic and molecular species 

A hydrogen atom impinging on graphene with a low 
collision energy can either physisorb or chemisorb. The 
physisorption regime has long been probed with the help 
of selective-adsorption resonances in H atom scattering 
off graphite^^. The extrapolated value for the physisorp- 
tion binding energy (~40 meV) to a single layer is in 
very good agreement with recent theoretical studieP^. 
Physisorbed species are highly mobile and easily desorb 
from the surface since they couple only weakly with the 
substrate. For this reason, chemisorption turns out to 
be more interesting for graphene electronic structure 
engineering. 

Chemisorption of single H atoms on graphite has 
been studied since the theoretical works of Jeloaica and 
Sidis^^ and Sha and Jackson^'^, who first showed that 
it indeed occurs if the substrate is allowed to relax. 
Among the four possible adsorption sites the hollow 
and bridge were found not binding while the two kinds 
of atop sites (with or without an carbon atom on the 
layer underneath for graphite) give essentially the same 
behaviour, since graphene layers in graphite lay ~ 3.4 
A apart. This implies that the (surface) chemistry of 
graphene is very similar to that of graphite. 
Adsorption on the top site induces a surface reconstruc- 
tion ('puckering'). Such a reconstruction consists in 
the outward motion of the carbon atom beneath the 
adsorbed hydrogen, and occurs as a consequence of 
sp^ — sp^ re-hybridization of the carbon valence orbitals 
needed to form the CH bond. The re-hybridization 
induces a change in geometry of the substrate site, 
from a planar {sp^) to a tetrahedral {sp^) form, thereby 
leading to the surface puckering. The energy required 
for such a process, defined as the energy difference 
between the relaxed and the puckered configuration, is 
substantial (~ 0.8 eV) and this explains why binding 
energies to graphene are typically smaller than for other 
carbon species. If the graphene layer is k ept fi at the 
carbon - hydrogen bond is metastable onl}^^'^, while 
allowing surface relaxation chemisorption becomes an 
activated process with stable products (~ 0.80 eV). 
When the hydrogen atom collides on a carbon site 
already puckered, i.e. already in the sp^ form, 
chemisorption is a barrierless process. Otherwise, 
following an adiabatic path (hence allowing the carbon 
atom relaxation to its equilibrium position at every point 
along the reaction coordinate) an energy barrier ~0.2 eV 
high is found, as a consequence of the re- hybridization. 
These adsorption curves are shown in figjS] The barrier 
(which is also present when the substrate is kept planar) 
has an important, purely electronic origin. Indeed, it has 
been shown^^ that it results from an avoided crossing 
between a repulsive interaction with the Kekule-like 
ground-state and an attractive interaction with the 
low- lying, Dewar-like excited state (see Fig. 4 in 130]) . 
This can be nicely understood in terms of the chemical 
picture above since the Kekule-like structures do not 
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FIG. 5. Left Panel: chemisorption curves for a H atom on a graphene top site as obtained from periodic DFT calculations 
using a 2x2 supercell. In the adiabatic path (black) the carbon atom is allowed to relax at each C-H distance, while in the 
diabatic path (blue) the top site is kept fixed in the puckered geometry, z is the adsorbate height on the surface. The inset 
on the bottom shows the equilibrium geometry. Right panel: map of potential energy surface as a function of the heights of H 
and C atoms above the surface, for a collinear geometry. Also indicated the diabatic and adiabatic paths. Units are eV and A. 



have unpaired electrons which can readily couple with 
that of the incoming H atom. 

The overall binding picture of H atoms has found sub- 
stantial experimental proof for graphite surfaces. Hot 
hydrogen atoms produced by dissociating H2 molecules 
at ^ 2000 K are required to overcome the barrier and 
observe sticking. Thus, chemisorption is under kinetic 
controi^^ ^ in marked contrast with vacancy formation 
through e~/ion bombardment discussed above. Indeed, 
as we show below, H atoms do not adsorb completely 
random on the surface. A number of TPD, AES, EELS 
and HREELS spectroscopy dats^^^ is available, along 
with detailed kinetic Monte Carlo simulations^^ of TP 
desorption curves and accurate studies of vibrational 
relaxation dynamics^^ and reaction dynamics to form 
^^40 46 Notice that even though we focused here on 
adsorption of H atoms the same holds for other simple, 
monovalent chemical species. 



IV. LOW DENSITY: tt-DEFECT STRUCTURE 

A. The appearance of midgap states 

The effect of atomic scale defects in graphite, and later 
on in graphene, has been experimentally studied since the 
late eighties, when scanning tunneling microscopy (STM) 
allowed to capture images on solid surfaces at atomic 
scale resolution. It appeared immediately that when a 
vacancy was created by irradiating the sample, a bright 
\/3xA/3i?30° charge density reconstruction appear d^^ . 
A carbon vacancy or a defect in the 7r-network due to a 
monovalent chemisorbed species creates in graphene an 
imbalance between the number of sites in each sublat- 
tices. This lowers the overall lattice symmetry, up to 



eventually remove the Dirac cones and open a band gap. 
Looking at the tight-binding Hamiltoninan in equation 
[1] the introduction of a 7r-defect in the graphene lattice 
reads as the removal of the basis function corresponding 
to the defect site, and the system eigenstates become 
necessarily odd-numbered. Therefore, in the nearest- 
neighbour approximation, because of the electron- hole 
symmetry, one of the eigenvalues in the energy spectrum 
necessarily lies at the Fermi level. This zero-energy state 
is a singly occupied molecular orbital called midgap state, 
even when a gap is not really present. When relaxing the 
nearest-neighbour approximation such state moves from 
the Fermi level, but remains close to it. Its presence is 
important for the transport properties being responsible 
for resonant scattering mechanisms. 
The appearance of midgap states in bipartitic systems 
has been intensively studied in solid state physics be- 
cause of the implications they have for the appearance 
of magnetism. Inui et al.^^ formulated a useful theo- 
rem for bipartitic tight-binding models with a sublat- 
tice imbalance. According to their result in any bipar- 
tite lattice in which the numbers of sublattices sites ua 
and ub are not equal, there are at least 77 = \nA — 
linearly independent eigenfunctions of the Hamiltonian 
at zero energy, all with null amplitudes on the minority 
sublattice sites. The proof is simple: for let Na > Nb 
and = Q^i \cii) be a trial solution at zero energy. 
The coefficients need to satisfy {bj\H\ai) ai = 
for j = 1^..Nb which is a set of Nb equations for the 
Na > Nb coefficients, with r] linearly independent solu- 
tions. This also shows that ip^s localize on the A lattice 
sites. 

Analogous results have been already known in hydrocar- 
bon chemistry for some time. The tight-binding approach 
described above has been used for decades in quantum 
chemistry to study aromatic hydrocarbons, under the 
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FIG. 6. Valence Bond model for the binding of a radical 
species (a H atom) on benzene. The numbers give the weights 
of the corresponding VB structure^^. The resulting unpaired 
electron localizes mostly in ortho (two leftmost structures on 
the top row) and para (mid panel) position, as emphasized 
by the bond-switching mechanism reported in the chemical 
formula of the bottom row. 



name of Hiickel method. The mathematical properties of 
the Hiickel Hamiltonian have been formalized in a series 
of theorems and corollaries in a famous book of Dewai^^. 
Bipartite lattices were listed there as "alternant" hydro- 
carbons, and the emergence of midgap states formally 
predicted in case of odd-numbered alternant hydrocar- 
bons. 

Calculations, both at tight-binding and at higher lev- 
els of theory (DFT), confirm these expectations: in 
graphene the zero-energy states originated in this way 
correspond to semilocalized modes around the defect 
which dec ay sl owly with the distance, i.e. with a 
power la W l ^^ l, a result which has been recently con- 
firmed by experiments^. Pereira et al.^^ performed a 
comprehensive analysis of the effect low-density defects 
have on the graphene DOS, by using numerical tight- 
binding calculations for ~ 4x10^ lattice sites and ana- 
lytic results. Analogous results have been found in DFT 
studies^ of isolated vacancies^^ and adatoms^^^*^. 



stitutions. In contrast to the full analysis of possible 
spin-couplings, exporting this model to graphene is then 
rather straightforward^. A picture of the mechanism is 
shown in figj7|for a coronene model, that is meant to rep- 
resent the whole graphene lattice. The itinerant electron 
hops between sites of one type only, thereby occupying 
a delocalized state which is the midgap state described 
previously in the tight-binding (MO) picture. 
Whatever picture we use the result is a spin density 
(magnetization) localized close to the defect, on the sites 
of the hexagonal sublattice not housing itP^. At low 
density, where hybridization does not occur, such spin- 
density thus determines the appearance of (microscopic) 
magnetically ordered domains. It further influences re- 
activity of the substrate with foreign species, which can 
readily 'saturate' (singlet-couple) this electron if they 
land on the correct sites, as will be shown in the next 
section. 

Before concluding this Section, we can now understand 
why simple adatoms do not move on the surface. Indeed, 
for the H atom to hop on the neighbouring site the un- 
paired spin has to move from one sublattice to the other 
and this requires breaking completely the existing CH 
bond and forming a new one: the barrier to diffusion, 
then, matches the desorption energy. This explains the 
experimental observation that H atoms are immobile on 
the surface^^. For more complex species, e.g. O atoms, 
spin-recoupling on the adatom may help the diffusion 
(isomerization) process via formation of a 'bridge' be- 
tween the two sites. This would explains why DFT com- 
puted barriers for diffusion of OH species are definitely 
smaller than the desorption energy^ 



V. HIGH DENSITY: SPIN-ORDERING, CLUSTERING 
AND RELATED ISSUES 



B. Chemical resonance formula 

In the case of a single 7r-defect, or a random distri- 
bution of them, the appearance of midgap states may 
be easily understood by applying the resonance-based 
VB picture described in Section |TT] . Considering ben- 
zene as the simplest building block of graphene, it is easy 
to realize how adsorption of a H atom breaks the aro- 
matic network and leaves one unpaired electron free to 
move on the lattice by bond switching: spin-recoupling 
with a neighbouring double bond creates an unpaired 
electron in one every two lattice sites. Ab-initio VB 
calculation^^ show that this indeed the case: the 5 tt 
electrons have 5 different ways of couplings (Figj6| but 
only those with the unpaired electron in the so-called or- 
tho and para positions are relevant; an electron in meta 
position would involve a Dewar-like structure, which has 
a high energy bond-pattern (see Figj6|. The bond switch- 
ing mechanism is very useful and well known in basic 
organic chemistry, where it easily allows predictions for 
orientation effects, e.g. in electrophilic aromatic sub- 



A. Predicting midgap states and magnetism 

The very simple counting rule for midgap states intro- 
duced above usually works fine for graphene, but fails 
to be predictive for some class of finite size graphenes or 
analogous (complementary) holes on the graphene sheet. 
For instance, the first two molecules in FigjS] have no 
sublattice imbalance (being symmetric) but are radical 
species, i.e. they necessarily have midgap states^^. To a 
closer inspection, sublattice imbalance is indeed only a 
sufficient condition for midgap states to appear. 
To rationalize the situation, it is necessary to introduce 
the concept of non- adjacent sites in a TV-site bipartitic 
system. We say that two sites are non- adjacent if they 
are not bound (connected) to each other; for instance, 
two sites on the same sublattice are non-adjacent. 
Clearly, there exists a maximal set of non-adjacent sites 
and we call a the sites in this set, and /3 the remaining 
ones {Noc^Nf3 = N — Na in number, respectively). 
Each site a binds at least to one site /3, otherwise it 
would represent a completely isolated site. Arranging 
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FIG. 7. Itinerant electron model for the -vacancy-induced midgap state 




FIG. 8. Molecules with no sublattice imbalance. The first two molecules on the left are di-radical species, i.e. they have two 
midgap states. The Schlenk-Brauns hydrocarbon shown in the middle panel differs from the Chichibabin hydrocarbon shown 
in the right panel by the connectivity only. Nevertheless, the latter does not present unpaired electrons. The largest set of 
non- adjacent sites is marked by red dots. 



one electron per site a, however, we can form at most 
N(3 bonds at a time, and therefore we are left with 
T] = Na —Nf3 = 2Na —N unpaired electrons, i.e. midgap 
states. The case of a sublattice imbalance discussed 
above is a special result of this rule: when Na > 
since the A lattice sites are always non- adjacent and 
^ we have r] > Na — Nb- Fig |8] shows molecules 
having > and no sublattice imbalance, with the 
indicated a sets. As it is evident from its derivation, 
this result can be equivalently re-phrased by defining 
77 to be the number of unpaired electrons in the Lewis 
structure(s) with the maximum number of tt {i.e. 
double) bonds. 

Notice that, since the spectrum of the Hamiltonian 
is determined by the system topology the whole set 
of counting rules for midgap states can be derived 
entirely from graph theory. In particular, midgap states 
appear as zeros of the characteristic polynomial of the 
adjacency matrix A, that defines the connectivity of the 
graplPIl^. In this context, the above result is known as 
graph nullity theorem. 

Having derived the exact conditions determining the 
appearance of midgap states, the question arises of how 



spins couple when a number r] of unpaired electrons are 
present. The determination of the spin state cannot 
come, of course, from the simple tight-binding Hamil- 
tonian, since in these open-shell configurations energy 
ordering is mainly determined by electron correlation. 
At first glance, it can be guessed that electrons oc- 
cupying quasi-degenerate midgap states tend to keep 
their spins parallel, in a sort of molecular Hund's rule, 
as this reduces Coulomb repulsion, i.e. system's total 
spin should always be^ rj/2. This is actually the case 
only when midgap states originate from a sublattice 
imbalance, since in such instance they are forced to stay 
on the same sublattice. When midgap states (unpaired 
electrons) lie on different sublattice they best couple 
at low spin. This result can be shown to be exact for 
the realistic model provided by the (repulsive) Hubbard 
Hamiltonian: LietP^ showed that for any bipartitic 
system at half-filling the ground-state spin S is given 
by the sublattice imbalance S = ^\Na — NbI- This is 
a subtle effect of electron correlation, which would lead 
to an energetically unfavourable spin polarization of the 
remaining occupied orbitals if the above Hund rule were 
followed in absence of sublattice imbalancd^. From a 
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different perspective, it has been associated with the 
most "spin-alternant" structure (Ovchinnikov's rul^^. 
According to the rules above it is now possible to predict 
the number of midgap states and the spin state of a 
number of complex graphene structures without relevant 
exceptions. We only note that the theorems stated 
above for bipartite lattices do not apply for topological 
defects that destroy bipartitism. Nevertheless, it has 
been noticed that the Ovchinnikov's rule can be usually 
extended to non-bipartite systems^^, although some 
care has to be paid^^. For instance, the ground- state 
multiplicity of Stone- Wales defects is correctly predicted 
to be zero by this rule. 



B. Preferential sticking 

When adsorbing hydrogen atoms on graphite or 
graphene under kinetic control STM images clearly show 
the formation of dimers and cluster^^. Since H atoms 
are immobile on the surface this must be due to a pref- 
erential sticking mechanism. This mechanism was first 
suggested by Hornakaer et al]^ who looked at the STM 
images formed by exposing Highly Oriented Pyrrolitic 
Graphite (HOPG) samples to a H atom beam, and ob- 
served formation of stable pairs, also confirmed by first- 
principles calculations^^ Later Casolo et al.^^ showed 
that the preference for certain lattice sites comes from 
the spin density localized on one of the two sublattices 
(the midgap state), as generated by the first adsorbate. 
The overall picture^^ is consistent with the VB chemi- 
cal model: when a first H atom is on an A-type site, 
the unpaired electron localizes on the B sublattice and 
bond formation easily occurs on its sites. An "AB dimer" 
(which has no sublattice imbalance) is formed and a sin- 
glet ground-state is obtained where aromaticity is par- 
tially restored. Conversely, if adsorbtion occurs on the 
same sublattice, i.e. to form 'A2" dimers, the incoming 
H atom does not make use of the available spin-density, 
and adsorption energies are comparable to that of the 
first H atom. Furthermore, as another electron is set free 
on the same B sublattice occupied by the unpaired elec- 
tron, the ground state is a triplet [r] = Nb — Na = 2). 
The results of DFT calculationP^on a number of dimers 
are shown in Figj9]as function of the site-integrated mag- 
netization, i.e. the average number of unpaired electrons 
in each site as results from the first adsorption event. It 
should be noted that substrate relaxation effects, though 
substantial (~ 0.8 eV), are site-independent for all but 
the ortho dimeJ^, thus the curves in the graphs of Fig. 
[9] refiects purely electronic effects. Binding and barrier 
energies both depend linearly on the local magnetization, 
thereby implying a linear relationship between them; this 
is a common tendency in activated chemical reactions 
known as Br0nsted-Evans-Polayni rule. An exception is 
provided by the ortho dimer (rightmost data point in the 
graphs of Fig. whose formation requires further rear- 



rangement in the first C-H neighbourhood. This is shown 
in figure [To] where the equilibrium geometry of the dimer 
is reported in the left panel. It is clear from the figure 
that the two H atoms point in opposite directions (as in a 
H-C-C-H eclipsed conformation of an alkane), which sug- 
gests that, despite their proximity, they would not easily 
desorb to form H2 upon heating the substrate. This is in- 
deed what has been found by a combined theoretical and 
experimental study by Hornalaer et al?^. upon heating, 
the ortho dimer prefers to isomerize to the para dimer, 
which dehydrogenates easier {i.e. at a lower tempera- 
ture). The highest temperature peak in the TPD spectra 
corresponds then to this isomerization process. The para 
dimer itself, whose equilibrium geometry is shown in the 
right panel of FigjlOj forms abundantly when exposing 
graphene to a H atom beam, since its formation is barri- 
erless (see figj9|. This forms the basis for the preferential 
sticking mechanism first suggested by Horenkaer et alW^^ 
which is here summarized with the results of fig(9j namely 
formation of AB dimers is both thermodynamically and 
kinetically favoured over formation of A2 dimers and sin- 
gle atom adsorption. 

Unfortunatel}^^^ , the preferential sticking mechanism 
above works only for dimers. Once an AB dimer is formed 
(A2 is unfavoured) there are no further unpaired elec- 
trons available, and no bias on the adsorption of addi- 
tional H atoms. This is confirmed by DFT calculatiorP^ 
on a number of larger A2B2, A2B, A3B1 and A3 clus- 
ters. As expected from the VB model above, adsorption 
of a third hydrogen atom to a stable AB dimer parallels 
that of the first H, with essentially no preference towards 
any specific sublattice position, and always produce dou- 
blet structures (M = 1/i^) (see fig 11). Similar conclu- 
sions hold when adding a third H atom to the (magnetic) 
meta dimer A2: adsorption on B lattice sites is strongly 
favoured {Ei^ind = 1-2 — 1.9 eV) and produces doublet 
structures (M = 1 /i^), whereas H atoms bind to A lat- 
tice sites with an energy ~ 0.7—0.8 eV and produce highly 
magnetic structures (M = 3 jub) (see fig 11). Energy 



barriers to adsorption follow the same trend: calcula- 
tions show that, with few exceptions, barriers to sticking 
a third H atom compare rather well with that for sin- 
gle H atom adsorption for the processes AB^A25 and 
A2 ^Aa, and may be considerably smaller for A2 ^A2B 
ones. These three-atom clusters, similarly to the single 
H atom, necessarily bias the adsorption of a fourth atom. 
The computed binding and barrier energies for this pro- 
cess have been found to compare rather well with the 
dimer values, actually they nicely fit to the same linear 
trends shown in Fig|9] Finally, all the considered A2B2, 
A2B, A3B1 and A3 clusters have been found to have 0, 
1, 2 and 3 unpaired electrons in their ground-state, re- 
spectively, in agreement with expectations {i.e. either 
the VB model or the Lieb theorem). 
Few exceptions to this picture are for compact clusters 
where substrate relaxation does play some role, see e.g. 
the structures on the right of fig. [TT] Compared to other 
trimers, these structures are favoured because of the sub- 
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FIG. 9. Left panel: Binding energies for secondary H adsorption as a function of the site- integrated magnetization (Msi), for 
AB ( circles) and A2 (squares) dimers. Both singlet (red) and triplet (blue) solutions are shown in red and blue respectively. 
Also shown the data point for single H adsorption (black diamond) and a linear fit to the data set (solid line). Right panel: 
corresponding barrier energies for secondary atom adsorption (ground-state only). Data point at M^/^O is for single H 
adsorption. 





FIG. 10. Equilibrium geometry for ortho (left) and para 
(right) dimer structures. 



strate softening occurring after formation of the para 
dimer, which is in a typical boat configuration (fig, 10). 
Such softening is expected to reduce the relaxation en- 
ergy needed for the binding of the additional atom, with 
a gain of some tenths of eV on the overall energetic 
balance^. This would explain why experiments at in- 
tense H atom flux do not find a random distribution of 
dimers, as would be expected on basis of electronic effects 
only, rather clusters made up of a number of atoms^^ 
Though a detailed analysis would require the knowledge 
of the adsorption barrier for a rather large number of 
clusters, the linear relationship shown above may help 
in making educated guesses on the basis of the binding 
energies only. 



C. A route to graphane? 

Graphane is a novel two-dimensional material, namely 
the fully hydrogenated graphene sheet. The structure is 
still bipartitic and each sublattice bears all hydrogens 
on the same side of the lattice plane (in meta to each 
other), in such a way to have a chair-like configuration 



as in cyclohexane. It is an insulating material with no 
TT electrons^^ that might have been recently obtained 
by simply exposing graphene to cold hydrogen plasms^. 
This result is very interesting in light of what has been 
shown in the previous sections. In order to produce 
graphane by simple hydrogen exposure it is necessary 
that either H atoms adsorb selectively on one sublattice 
only for a given graphene face or hydrogen diffuses 
(isomerize) to occupy the sites on the right face. 
According to the discussion of previous section, forma- 
tion of meta dimers is unlikely, ortho, para positions 
being highly favoured for adsorption. According to DFT 
calculations on free-standing graphene^^, this is true 
both for the syn- (on the same face) and the anti- (on 
opposite faces) dimers. However, even if hydrogen atoms 
were likely to form anti-para dimers also in supported 
graphene (and they are required for graphane produc- 
tion) an efficient syn- to anti- conversion mechanism 
would be needed to convert those dimers already formed 
on the same graphene face. Unfortunately, as we have 
seen in the previous section, even in this fortunate 
case, no true preferential sticking can occur after dimer 
formation and it will be very unlikely that all the other 
hydrogens will chemisorbed in the correct sites and face. 
Indeed, recent molecular dynamics simulations showed 
that disordered, frustrated hydrogenated domains would 
rather form^^. 

If graphane has been really formed by hydrogen exposure 
some other effect has to play a role. Curved graphene 
areas might help this process. Graphene is a very elastic 
membrane that naturally exhibit ripples that tend to lay 
down along steps and kinks of the supporting substrate 
on which it was grown. Indeed recent experimental 
findings suggested that hydrogen chemisorbs more 
efficiently on the ridges of the silicon carbide substrate 
surface onto which graphene usually lays^A Moreover, it 
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0.820 eV 

FIG. 11. Some of the possible hydrogen trimers structures obtainable from the AB para dimer (left diagram) and from the A2 
meta dimer (right). Binding energies are also shown. 



has been also shown that hydrogenation of single-layer 
graphene is easier than for many-layer graphene, likely 
as a consequence of t he hig her corrugation displayed 
by the graphene surfac^^I'^. This is reasonable, as for 
nanotubes the curvature reduces the Pz-Pz overlap, i.e. 
aromaticity, thereby lowering the barrier energies for 
H chemisorptioil^^^. Still, there is no clear evidence 
that local curvature plays a role in graphane forma- 
tion, and more investigations in this direction are needed. 



VI. DEFECT-BASED MATERIAL DESIGN 

When it comes to device fabrication only few of the 
many extraordinary properties of graphene are relevant, 
at least for the chip-makerJ^. Among them, its thick- 
ness allows the thinnest possible gate-controlled regions 
in transistors and, according to scaling theories, should 
reduce electrostatic problems if short channels have to be 
built. Mobility is an important factor as it allows for in- 
stance high-performance interconnects and fast response 
to external (gate) potentials. It becomes of secondary 
importance in short channels, where high fields build up 
and carrier velocity saturates, but also in this respect 
graphene proved to have superior properties than con- 
ventional materials. Indeed, high-performance transis- 
tors for frequency applications have been realizecP^, and 
record cut-off frequencies are being continuously scored. 
However, for its usage in logic applications the absence of 
a band-gap is a major probler d^^^ even when the Fermi 
level crosses the charge neutrality point a non-zero resid- 
ual conductivity avoids the complete current pinch-ofl"^^ . 
The absence of a band-gap, indeed, prevents the achieve- 
ment of the high current on-off ratios required for logic 



operations. 

Graphene can be turned into a true semiconductor by 
properly engineering it. Electron confinement, though in 
general not trivial for massless, pseudorelativistic carri- 
ers, can be obtained by cutting large-area graphene to 
form narrow nanoribbons. Apart from related fabrica- 
tion issues, one main drawback of such an approach is 
the removal of the Dirac cones and the resulting band- 
bending. This is expected to increase the effective mass 
of the carriers, thereby reducing their mobility. Indeed, 
it has been generally found that mobility is a decreasing 
function of the gap^^, and this is an undesirable by-side 
effect worth considering with such a traditional approach. 
Alternatively, symmetry breaking is known to turn the 
massless Dirac carriers into massive (yet pseudorelativis- 
tic) carriers. This can be realized by depositing or grow- 
ing graphene on a substrate that renders inequivalent 
the two sublattice positions. For instance, boron nitride 
has the same honeycomb lattice as graphene and a simi- 
lar cell parameter, but presents two inequivalent sublat- 
tices. When graphene is in contact with such a surface B 
and N interact differently with the carbon atoms of the 
graphene sheet, breaking its sublattice equivalence and 
lifting the degeneracy of the two bands. A similar situ- 
ation is achieved for graphene grown on silicon carbide 
surfaces, where a gap has been observed by angle-resolved 
photoemission spectroscopy though subtle electron cor- 
relation effects may play a role in such cas d^^ * ^^ l 
In the following sections we describe alternative possibili- 
ties for opening a gap in graphene band structure, namely 
those offered by superlattices of defects and dopants. 
One interesting finding in this context is the proof that 
a band-gap can be opened in graphene without break- 
ing its symmetry, with the advantage the new Dirac 
cones (massless carriers) appear right close to the gapped 
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regioiPI'. For this reason we start introducing some sym- 
metry considerations, extending the arguments given in 
section [III 



A. Symmetry considerations 

As we have seen in section [IT] graphene's unconven- 
tional electronic properties are strictly related to its 
Dq^ point symmetry. The A:-group at the K-K' high- 
symmetry points (Dsh) allows for doubly degenerate ir- 
reducible representations, and Bloch functions built with 
Pz orbitals of A and B sublattices span just one of its 
two-dimensional irreps. As e — h symmetry does not mix 
one- {A) and two- (E) dimensional irreps this level has to 
lie at zero energy, where the Fermi level {ep) is located. 
Were not there such degenerate level, graphene would 
be, as any other bipartitic system at half- filling, semicon- 
ducting. Graphene can be forced to be so by either low- 
ering the symmetry (i.e. changing the k group at K(K') 
to a simpler one), or changing the number of E irreps 
at the special points while keeping the overall symme- 
try. In the latter, more intriguing case, since the overall 
point symmetry is preserved, degeneracies may still oc- 
cur at energies different from e^?, and new Dirac cones 
are to be expected. The "recipe" for doing that is very 
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simple^ nxn graphene superlattices have the same sym- 
metry properties and 2n^ atoms per cell; by symmetrical 
removing a number of C atoms^^ is possible to change the 
number of irreps and turn, in particular, the E ones to be 
even at every, highly symmetric special point (T^K^K'). 
With few exceptions of residual accidental degeneracies, 
this opens a gap in the band structure. 
The approach is made effective by counting the number 
of irreps generated by 2n^ atoms in a nxn unit super- 
cell. The results of this calculatiorP^ can be grouped 
into three different sequences, n = 3m, 3m + 1,3m + 2 
(m integer), according to the BZ folding properties. In 
two thirds of the cases, i.e. when n = 3m + 1, 3m + 2 , re- 
moval of the atoms at the center of the two-half cells (red 



balls in the left panel of fig 12) is sufficient for opening 
a gap. Figure 12 (right panel) shows one of the simplest 
resulting semiconducting superlattice, namely with sim- 
ple, atomic-scale defects arranged in a honeycomb lattice. 
Because of their nature, they are best considered as su- 
pergraphenes. 



pare favourably with the gap in armchair nanoribbons^. 
However, one distinctive feature of such structures is the 
additional presence of new Dirac cones right close to the 
gapped region (blue circle in figjl3|. This might be im- 
portant in charge transport, since they can sustain mass- 
less carriers when the Fermi level, as tuned by a gate 
potential, is swept across the gap. 

In practice, it is still experimentally challenging to 
realize the atomic-scale patterned structures introduced 
above. It is however sufficient to consider similar 
superlattices of holes analogously to the graphene 
antidots superlattice j^^^^ . The resulting structures are 
honeycombs antidots as the one shown in figjM] They 
are experimentally feasible, since it has been shown 
possibl d^^tiiQi l to create circular holes with diameters as 
small as 2 — 3 nm and periodicity 5 nm. Analogous 
patterns of H adatoms have also been realized thanks 
to the interaction between graphene and an underlying 
metal surface that creates Moire patterns activating 
chemisorption in specific area^^^. 

Tight-binding calculations on the honeycomb antidots 
show that the band-gap is quite large for reasonable 
values of the superlattice constant and of the hole 
diameter and, as before, new Dirac cones appear at low 
energies, close to the gapped region^. The gap size is 
even larger than for the simplest structures considered 
above, though the latter remain optimal in this context: 
when the gap size is renormalized to the number of 
defects per unit cell, honeycombs such as that reported 
in fig. [T2] prove to be "magic". 

Notice that previously suggested antidot 
superlattices^^ show comparable gaps, and are 
therefore equally valid candidates for turning graphene 
into a true semiconductor. The only difference is a subtle 
symmetry-related issue. These hexagonal superlattices 
are all of >/3nx>/3n type, hence with the same Dq^ 
symmetry of the honeycomb lattices considered here, 
and this would suggest that exactly the same results 
hold for them. A closer inspection, however, reveals 
that in the >/3nx>/3n case, the K, K' points of pristine 
graphene always fold to F. This is advantageous for the 
band-gap opening, since these structures are generally 
semiconducting if sufficiently defective^^^. However, 
as in r the k group has inversion symmetry, residual 
degeneracies at e 7^ e^?, is not lifted at first order, and 
therefore no linear dispersion is present. 



B. Superlattices of vacancies or holes 

Tight-binding and DFT calculations on the simple 
structures identified in the previous section show indeed 
a sizable band-gap. The gap size approximately scales as 
vp/ln where is the Fermi velocity in pristine graphene 
and In is the distance between defects {In = na/VS), 
as can be guessed from a dimensional analysis or ob- 
tained from a perturbative calculation within the tight- 
binding approach^^. Both the size and the scaling com- 



C. Superlattices of substitutional atoms 

Symmetry arguments similar to the one given above 
apply as well to the case where C atoms are replaced 
(rather than removed) by other species, in such a way to 
form superlattices substitutional dopants. The only dif- 
ference is that now foreign species are present and point 
symmetry can be altered. Here we focus on group IIIA 
and VA elements, mainly because of the fast progresses in 
methods for the controlled synthesis of N- and B- doped 
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FIG. 12. Left: counting the number of irreps generated by the atomic basis in a 4x4 supercell: indicated are the irreps 
generated by the atoms at the center of the half-cells (red balls) and by green triangles. Right: a simple "supergraphene", the 
simplest defective 14x14 honeycomb. 
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FIG. 13. Energy gaps in simple supergraphenes made with pz vacancies. Left: results of TB calculations as functions of 1/n. 
The symbols are for different parameters of the TB hamiltonian, and the solid line is the result of a perturbative calculation 
at the K point. Se^^ for details. Middle: tight-binding (red) ^5. DFT (black) results. In the latter case, defects have been 
modelled as H atoms. Stars represent the results for asymmetric dimers placed in the same nxn supercells. Right: energy 
bands for the n = 13, 14 supergraphenes. 



graphenes. For instance, Panchakarla et aP^ have re- 
cently shown how it is possible to insert B or N dopants 
in graphene by adding the correct precursors in the arc 
discharge chamber, while Ci et al.^^^ have reported the 
synthesis of large islands of boron nitride embedded in 
graphene by atomic layer deposition techniques. Meth- 
ods to selectively replace C atoms from graphene lattice 
have also been proposed by Pontes et al^^^. 
Substitutional defects behave similarly to Pz vacancies 
(to which they reduce when the hoppings become zero) 
but introduce impurity bands which partially hybridize 
with those of the substrate. In addition, the diagonal 
disorder they introduce breaks e — h symmetry giving 
rise to a Fermi level shift, i.e. to p— and n— doping for 
group III A and VA elements, respectively. If superstruc- 
tures are only weakly defective, however, the Fermi level 
shift scales as 1/n, since the linear-energy dispersion im- 



plies Ep = VF\/^Tie (here rie is the electron (hole) excess 
density, ng oc l/v?). Thus, analogously to the superlat- 
tices of the previous section, the defect-induced pertur- 
bation affects the electronic structure close to the Fermi 
level, and symmetry arguments may be used to establish 
whether degeneracy occurs at the special points in the 
important low-energy region. 

It has been shown that, depending on the overall sym- 
metry, superlattices of N and B defects (and mixture 
thereof) can either preserve the Dirac cones {Dq^ super- 
lattices) or open a band gap {D^hV^. For instance, hon- 
eycomb superlattices of B (or N) dopants only (I^e/i sym- 
metry) are found to preserve the Dirac cones, and to be 
only p— (n— ) doped (the cone-apex shift being oc 1/n). 
Indeed, the low-energy band structure in these superlat- 
tice is very similar to that of graphene, with a group 
velocity at the shifted Dirac cones depending linearly on 
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FIG. 14. A honeycomb antidot superlattice presenting a gap 
by the symmetry preserving approach discussed in the text. 



1/n too. As the Fermi level can be tuned by a gate po- 
tential, these systems offer the possibility of investigating 
the role that the effective speed of light (the above veloc- 
ity at the Dirac cones) has on the transport and optical 
properties of graphene. Conversely, BN-honeycomb su- 
perlattices {Dsh symmetry) , thanks to the balanced dop- 
ing, are found to develop a gap right at the Fermi level. 
The resulting dispersion relation is found to be quasi- 
conical, corresponding to massive Dirac fermions. The 
resulting gaps are found sizable and, furthermore, the 
effective rest masses are rather small, m* ~ O.Olmg or 
smaller. This suggests that these structures might have 
good field-switching and transport properties. 
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